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Nutation Damper for the AMPTE-IFM Satellit e 
Final Report 


SUMMARY 

This report describes the design of, computations on and testing of the 
nutation damper for the AMPTE-IRM satellite. A passive liquid damper has been 
developed building on the investigations carried out during the feasibility 
study. It consists of a liquid filled O-shaped tube system. The nutational 
motions induce liquid oscillations in the tubes; energy is destroyed by 
tube friction, constriction of the stream and the formation of turbulence 
which naturally is extracted from the satellite movement which is thereby 
damped. 

The equations of motion (non-linear coupled differential equations 
for satellite rotation velocities and the deflections of damping fluid) are 
derived for determination of the design parameters; using simulators the 
oscillation behaviour can be assessed and the parameters fixed. This design 
represents a good compromise with respect to the twelve different satellite 
configurations. The attenuation behaviour is evaluated using time constants 
( reduction in amplitude to 1/e). 

Tests carried out in a pendulum testing device to check out the theory 
and to determine friction loss numbers show the anticipated dependency of 
the time constants on the respective excitation (= nutation) frequency 
and good conformity with the theory. 

Along with this, the time constants for the satellite are calculated. 

For the design selected (overall length = 572mm, weight = 1165 g, mounting 
radius = 375mm) values of 5.6 to 12.3 min result; they are significantly 
under the specified value of 30 min. Therefore, from the standpoint of the 
nutational motions, the reliable execution of the mission is assured. 


This study was prepared by Dr. A. Truckenbrodt (Project Manager) , Dr. B. Scholt- 
ysik and Dr. J.P. Mehltretter (scientific advisor) under the overall respon- 
sibility of J.M. Mehltretter. 

* Numbers in the margin indicate pagination of foreign text. 
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1. INTRODUCTION 


tv 


The intended objective was the design, manufacture and testing of a 
nutation damper for the AMFTE-IRM satellite. In accordance with the results 
achieved in the feasibility study, a passive fluid damper was developed. 

The emphases of the tasks described in this final report are 

a) adaptation of the various damper parameters ( dimensions, fluid, etc) to the 
specified conditions, that is to: 

a 12 different configurations of the satellite ( with respect to moments 
of inertia, rotational velocity), 

• maximum permissible damping time constants; 

b) experimental testing of the theoretical data and the determination of 
the import an ' damper parameters. 

The satellite configurations and the most important requirements, ref er- 
ed to the feasibility study, are brought together in the following: 
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maximum nutation angles: 


+ 20° for configurations 1,2 
+ 3° for configurations 3-12 


residual nutation angle: 
time constants; 
weight 


+ 0 . 2 ° 

r < 30 minutes, 15 minutes preferred 
k 1 kg/damper 


2. DESCRIPTION OF THE NUTATION DAMPER 


15 / 


2.1 Functional Principle 

The nutation damper shown in Figure 1 consists of a closed tube system 
that is constructed of two long thin tubes ( "tube") and two short thick 
tubes ("cap"). The damper is partially filled with a fluid. Based on the 
centrifugal acceleration through the rotation ( ) around the satellite spin 

axis, the fluid is located in the tube the furthest removed from the spin 
axis and in the adjacent areas of the caps. The nutation motions of the 
satellite excite fluid oscillations in the damper tubes; kinetic energy is 
destroyed (i.e. converted to heat) through the following effects: 

• tube friction of the flowing liquid in the tube 

• constriction of the stream by cross-sectional change from 
cap to tube 

• formation of turbulence by stream enlargement through cross- 
sectional change from tube to cap. 

This energy loss is extracted from the satellite such that the nutation 
motion is reduced. 



Figure 1 Passive Fluid Citation Damper 
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The effectiveness of the damper increases with the amplitude /6/ 
of the fluid oscillation. For this reason, the best damping effect is achieved in 
the case of resonance tuning (i.e. the nutation frequency v = characteristic 
frequency of the fluid oscillation ) , see Figure 2 . Moreover, one can 
see from Figure 2 that in case of increasing damping the effectiveness of the 
nutation damper 

• falls off in proximity to the resonance point ( v/ui p * 1 , corresponds 
to configurations 3-12) 

• increases for frequency ratios v /^ « i (corresponds to 

F 

configurations 1 and 2) 

• is less strongly influenced by the frequency ratio V//uj p . 

2.2 Damper Design and Structural Configuration 

In order to approach resonance tuning as closely as possible for configu- 
rations 3-12, a relationship for the fluid characteristic frequency to 
the spin frequency w zo of 


t V a) zo = °* 4928 

was used. This frequency relationship, given a mounting radius r Q = 375nm 
is achieved with the following damper geometry: /7 / 

Tube length = 472mm 

Tube interior 

diameter = 19nm 

Cap length = 80mn 

Cap interior 

diameter * 48.6mm 

Figure 3 shows the structural form of the nutation damper. Since a 
non-magnetic material is required , and on weight grounds as well, the 
damper was made of aluninun. Off-the-shelf aluninun tubes were used for the 
tubes. Due to the low wall thickness, the tubes and caps were not Joined by 
welding but by adhesives. For this reason, the caps are fabricated complete 
with the couplings needed for the adhesion. 

The filling of the danper with the fluid takes place at ambient pres- 
sure so that in space an overpressure of 1 bar prevails. Before gluing on 
the cover, heliun is piped into the damper so that a hellun tightness test 
can be carried out. An aluninun fixture serves to anchor the damper in 
the satellite. 
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Figure 2 : Influence of the Frequency Ratio v/u F and 
of the Damping on the Effectiveness of the 
Nutation Damper 
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In the case of the nutation damper design (Figure 3): 

• Overall weight (filled damper with mounting) 

• Overall length 

e Cap outside dimensions: Diameter 

Length 


1165 g 
572.2 ran 
51.6 mn 
90 nro 


The fluid "Flutec PP3" is used to fill the damper: 

• Fill level in cap: « 40 mn 

•» 

• Fluid volume: » 282.2 cm 




: Structural Configuration of the Nutation 
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3. THEORY 


N 


3.1 Equations of Motion 

The equations of motion for satellite and damper are contained In the 
feasibility study In linearized form. For an accurate estimate of the Influence 
of the non-linearities, the general equations of motion were derived for the 
coordinates 



(1) 


(b> - satellite rotational velocities, z^, fluid deflection In tubes 1,2); 
In general written form they read: 


M <*) i 


+ L 

'll 


7 .) 


mit E 


0 0 0 
0 0 0 


- 0 



( 2 ) 


(3) 


Here, in addition to the laminar tube friction, energy losses 
were also taken into accowt vrfiich are caused by the in and out streaming from 
tube into cap (stream constriction, turbulence formation) and turbulent tube 
friction (alternatively to the laminar tube friction), all in the interest of 
broadening the description of the energy removal used in the feasibility study. 

For the energy loss of a damper it is thus found, see equation 
(13) in the feasibility study, that 


dr; 

cit" 


a t * 





(4) 


with the loss term z - , vrtilch is only experimentally determinable. The damper 
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equations must therefore, in accordance with equation (13) of the feasibility 
study be supplemented by the terms 


. t z* • Rgn 


Estimations with C « 2 show that the energy loss characterized by t is /: 
excessive at large nutation amplitudes. It decreases quadratically with 
decreasing amplitude in favor of the linear term of the laminar tube friction. 
3.2 Time Constants 

On account of the non-linearity of the equations of motion, the amplitude 
history is no longer characterized by an exponential function with constant 
exponents. Therefore, it is impossible to define time constants independent of 
the amplitude. Equally, a sufficiently accurate energy-sink approximation is 
no longer possible. Therefore, the "theoretical" time constants are ascertained 
by the following procedure: 

1. Simulation of complete non-linear equations of motion 

2. Conversion of the amplitude term A A after an "examination period" 

At into an equivalent exponential function or time constant 

AA « e " At/x 

3. Time constant 

, T *= - At/ In AA 

The time constants derived with this procedure are collated in Section 5. 

4. TEST DEVICE 

Damping tests were carried out with 

• plexiglass nutation dampers with various filling levels 

• the ultimate flight model 

to check the theory and to experimentally determine the loss term C ; these 
were done using the test device described in the feasibility study. With the 
plexiglass dampers, the stream processes could be clearly observed. 

The test mechanism ( Figure 4 ) consists of a pendulun (total mass 
suspended at point A. By displacement of the pendulun tube (m^) and of the 
weights ny,, the center of gravity interval s and the moment of inertia 
0 A can be changed with respect to the support point A. In this way the period 
of oscillation T of the pendulun can be changed. 
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Figure 4 : Test Mechanism for Nutation Damper 

The nutation damper ND is fixed on the pendulum at interval a from the 
pivot A. 

At different pendulun geometries, (i.e. different characteristic frequencies 
of the pendulun), the oscillation behaviour for an initial amplitude of 
2 0 = 15° was recorded with a rotary potentiometer and recorder both 

with and without the nutation damper. Without the nutation danper the 
effect of bearing friction and air resistance is determined. The time constants 
amount here to about 200 to 850 sec depending on the pendulun moment of 
inertia. In the measurements with the nutation damper, the interval a 
was set at: a = 375mm = r Q , corresponding to the installation in the 
satellite. 

From the time constants measured for air resistance and bearing friction /12/ 
t (in the case of constant pendulun geanetry) and the time constants 
measured for the nutation danper + air resistance + bearing friction 
t follow the time constants T p of the pendulun only with 

the nutation danper: 


t nd+l+l * t l+l 
t l+l ~ t nd+l+l 


( 8 ) 


Figure 5 shows the time constants of the pendulun thus determined, rel- 
ative to the potential initial energy of the pendulun. 
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( 9 ) 


It is recognized that the effectiveness of the nutation damper is at 
its greatest (i.e. smallest time constant) when the pendulun (= excitation) 
frequency f^ coincides with the characteristic frequency fg of the fluid 
oscillation. The characteristic frequency on earth amounts to 




4 IT V 


FL 



2 


_) 

T 


f E = 23,76 oscillation/min 


In contrast, in space there follow: 

in configurations 1 and 2 : f E = 29.568 oscillations/min 

in configurations 3-12 : f^ = 9.856 " " 

Thus, in space the relationship between the nutation frequency v and /13/ 
the characteristic frequency f £ amounts to 

in configurations 1 and 2 : v /f^ = 0.41-0.48 
in configurations 3-12 : v /fg = 0.89-1.09 

These frequency ratios w re achieved on the test mechanism (on earth) 
with the pendulun frequencies 

fp = 9.7- 11.4 osc./min. £ configurations 1 and 2 

fp = 21.1 - 25.9 osc./min. £ configurations 3-12 
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pendulum frequency V 

•* • - Wk 

o.? c /. o r- c.e. c i r, 1.2 


corresponds to 

configuration 1,2 configurations 3-12 
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Figure 5 : Time Constants of the Test Mechanism Pendulun, relative to the 
potential initial energy, as a function of pendulun frequency. 
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Both the identically constructed nutation dampers of the satellite ("ND 1" /14/ 

and "ND 2" were used in the measurements on the test mechanism. It is recognized 
that the findings on both nutation dampers correspond well, as would be 
expected. Moreover, the limited deviations indicate that measurement and 
evaluation errors in the investigations on the test mechanism are also small. 

To determine the loss term L , the experimentally determined time constants 
for the test mechanism were conpared with the time constants computed for 
the test mechanism. Here it is shown that with L =2.2 a very good agreement 
is achieved between computation and measurement ( error about 1.5-4% ) in the 
neighborhood of the resonance point, that is for fp/f g = 0. 9-1.1. For the 
frequency ratio f p /f E = 0. 4-0-5, an optimun value L = 4.9 follows. These 
different L values make it clear that in the resonance case and with 
subcritical excitation, different stream conditions are present in the nuta- 
tion damper. 

5. TIME CONSTANTS FOR THE SATELLITE 

With the nunerical values that apply to the satellite, the time constants 
shown in Figure 6 result, as a function of satellite configuration. An initial 
nutation angle of =20° was used in configurations 1 and 2 for the com- 
putations and of ? f> = 3° in the succeeding configurations 3-12. The values 
ascertained experimentally on the test mechanism were used for the loss 
term t , that is: 

configurations 1 and 2: C =4.9 

configurations 3-12 : L =2.2 



| A > ■ — — — 

w^Orpm u» zo =20rpm 

A<T 20 ’ V 3 ’ 


Figure 6 : Computed Time Constants for Different Satellite Configurations 
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It Is seen that in all configurations the time constant lies between 
about 5 and 12 minutes. In configurations 1 and 2, the higher C value 
produces an improvement, because at a subcritical excitation the effectiveness 
of the nutation damper is improved with an increase in C (see Figure 2 ). 

For conpari >on, computations were also carried out in configurations 1 and 2 
using C 2.2. In this way the range of time constants vrtvLch result during 
the reduction of a large initial nutation angle down to a small final angle 
is limited. For configurations 3-12 (i.e. for v/f E ~ 1) one can expect an 
even faster damping at very small nutation angles because then the influence of 
the turbul* ice damping is smaller and near the resonance point a more effective 
damping is ichieved with the tube friction (see Figure 2 ). 

In Fi gure 7 are shown the final nutation angles after a total damping 

time of t = 1000 sec, vriien the initial nutation angle in configurations 
8 " 0 _ 

1 and 2 $ = 20 and in configurations 3-12, 6 q = 3°. 

/16 / 



1 2 3 4 5 6 7 6 9 10 11 1 ? 


Satellite Configuration 

Figure 7 : Nutation Angle of the Satellite after 1000 sec for the Different 
satellite Configurations 
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Figure 8 shows for the different configurations the maxinun initial 
nutation angle vd\ich results in a final nutation angle of B.^goo sec = o 2° 

after a damping time of 1000 sec. 



Figure 8 : Maximum Initial Nutation Angle with which after 1000 sec a 

Nutation Angle of 0.2° is Achieved, for the Different Satellite 
Configurations 
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Figure 9 shows , using the relatively unfavorable example of configuration /17/ 
5 (see Figure 6 ). that with a change in the characteristic value C , the com- 
puted time constant of the nutational motion of the satellite remains about 
in the same order of magnitude. This means that the time constants calculated for 
the satellite remain approximately correct, even when the actual characteristic 
value £ occurring in space deviates somewhat from the value developed on the 
test device. 



0 L V 

0 1 2 3 C 

Loss term £ 


Figure 9 : Influence of the Loss Term £ on the Time Constant of the 
Satellite (Example in Configuration 5) 
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With this it is shown that the damper clearly remains within the originally 
specified values, also estimated in the feasibility study, and accordingly 
it can provide for good dynamic behavior of the satellite. 
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